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In this review we discuss the current status of research in nuclear physics which is being carried 
out in different centers in the World. For this purpose we supply a short account of the development 
in the area which evolved over the last 9 decades, since the discovery of the neutron. The evolution 
of the physics of the atomic nucleus went through many stages as more data become available. We 
briefly discuss models introduced to discern the physics behind the experimental discoveries, such 
as the shell model, the collective model, the statistical model, the interacting boson model, etc., 
some of these models may be seemingly in conflict with each other, but this was shown to be only 
apparent. 

The richness of the ideas and abundance of theoretical models attests to the important fact that 
the nucleus is a really singular system in the sense that it evolves from two-body bound states 
such as the deuteron, to few-body bound states, such as “^He, ^Li, ®Be etc. and up the ladder to 
heavier bound nuclei containing up to more than 200 nucleons. Clearly statistical mechanics, usually 
employed in systems with very large number of particles, would seemingly not work for such finite 
systems as the nuclei, neither do other theories which are applicable to condensed matter. The 
richness of nuclear physics stems from these restrictions. New theories and models are presently 
being developed. Theories of the structure and reactions of neutron-rich and proton-rich nuclei, 
called exotic nuclei, halo nuclei, or Borromean nuclei, deal with the wealth of experimental data 
that became available in the last 35 years. Further, nuclear astrophysics and stellar and Big Bang 
nucleosynthesis have become a more mature subject. Due to limited space, this review only covers 
a few selected topics, mainly those with which the authors have worked on. Our aimed potential 
readers of this review are nuclear physicists and physicists in other areas, as well as graduate students 
interested in pursuing a career in nuclear physics. 
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I. NUCLEONS, NUCLEAR STRUCTURE, AND POLITICS 

Nuclear Physics evolved in a rather logical mannen from simple ordered shell model tailored in accordance with 
atomic structure to reproduce the magic numbers [l|, Q , to collective models needed to describe vibrations and cases 
involving non-spherical nuclei. The collective model is now considered an elaborate extension of the shell model 
through the inclusion of nucleon-nucleon correlations, which are predominantly of the pairing plus quadrupole type 
Q. In fact the shorter-range pairing correlations were introduced into nuclear physics at almost the same time that 
the BCS theory of superconductivity was developed [^. The quadrupole-quadrupole interaction is responsible for 
long-range correlations and gives rise to deformed mean field necessary in the description of deformed nuclei with 
their vibrational and rotational spectra Q. For general texts on nuclear physics theory we suggest Refs. [1, [^-[l^ 

Correlations involving particles and holes are also used in the description of excited states in nuclei, the most notable 
of these is the so-called giant resonances^ similar to plasmon resonances in atom clusters [IM3. In the 1970’s algebraic 
models based on group theory were introduced, guided by the pairing plus quadrupole models alluded to above. The 
resulting Interacting Boson Model (IBM) has as building blocks pairs of correlated particles (neutrons or protons) 
coupled to angular momentum 0 and 2 and treated as bosons [3, 1^ . . The resulting Hamiltonian can be analyzed 

using symmetry concepts and the corresponding spectrum is obtained analytically with the help of group theory. 

With the development of models of nuclear structure, which pointed out certain degrees of freedom to be studied 
to test these models, there was the accompanying development of a theory of nuclear reactions, aimed to probe these 
different degrees of freedoms. The concepts of compound nucleus [1, [23 - 1^ and direct reactions [l^ became a 
common language with the bombarding energy of the probing projectile being the deciding factor in applying either 
one of the two concepts: low energy projectiles are captured by the target to form a highly excited compound nucleus 
which then decays into the different open channels, while higher energy projectiles excites the target nucleus or transfer 
a nucleon or more Ira. A unifying concept that deals with both componnd and direct processes in an average way is 
the optical model |23l| , where the scattering Schrodinger equation is solved with a complex potential whose imaginary 
part is meant to simulate processes that remove flux from the elastic channel described by this same equation. It 
is this rich variety of possibilities which prompted the great advances made in the quantum scattering theory that 
underlines nuclear reactions. In this review we give a brief account of these development. 

As decades passed by, the main interests in nuclear physics were broadened and strong overlaps with other areas 
of science emerged. Our knowledge about the nucleons also changed with time. Now we know that they are built by 
deceiving confined particles] the quarks and gluons, generically known as partons. New theories and experiments are 
now dedicated to probe the sub-nucleon degrees of freedom and how they change in the nuclear medium, in stars and 
in the evolution of our Universe. We are witnessing an increasing enthusiasm in the nuclear physics community with 
its stronger involvement in problems with interface with many other areas, such as Chemistry, Biology, Cosmology, 
Astrophysics, Energy, Medicine, etc. 

Needless to mention that Nuclear Physics carries a big burden of being tightly involved with political issues such 
as national security and non-proliferation of nuclear weapons, a serious threat to the survival of the human race. 
But, it is also rewarding to witness the involvement of many nuclear physicists in politics, a good example being our 
colleague E. Moniz, a nuclear theorist and since 2013 the U.S. Secretary of Energy. Moniz was pivotal in worldwide 
non-proliferation negotiations. He has been regularly in Brazil during the 1980’s and 1990’s, e.g., during meetings 
organized by the present authors. He contributed to numerous important theoretical works in Nuclear Physics [13, [Ig. 

The stories we tell in this brief review are divided into two main topics: the nuclear physics proper, and nuclear 
astrophysics. Everyone says that figures are worth a thousand words. But due to the limited space, we will only use 
words (and some equations, as we are physicists, OK?) to tell our stories and point to a few open questions in this 
intriguing, beautiful, and extremely difficult field of science. 

In Section H we present an account of the different facets of nuclei and discuss the major advances in nuclear 
structure and put special emphasis on the recent research activities in the area of neutron- and/or proton-rich nuclei, 
which have been produced and studied since the 80’s. In Section HI we discuss reaction studies and discuss the 
different theories developed for the purpose. In Section IV we discuss the new nuclei we have discovered in the last 
few decades. In Section V we discuss the importance of nuclei in the Cosmos, and in the synthesis of elements. Einally, 
in Section VI, concluding remarks are given. 


II. DIFFERENT FACETS OF NUCLEI 

A. Quantum Mechanics, Nuclei, and the Origins of Life 

Nuclear Physics is the driving science behind stellar evolution. It determines energy production, element abundance 
and has ultimately led to the knowledge development in other basic sciences such as Astronomy. Imagine a world 
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without Astronomy and the benefits its observations have brought to such mundane things as the human sense of 
time, distance and motion, through navigation, and the discoveries of new worlds on Earth and in the skies. The 
symbiosis between nuclear physics and astronomical observations has skyrocketed in the last decades with the launch 
of astronomy dedicated satellites and the development of new observation technology. In the end, Nuclear Physics is 
the crucial piece of science to interpret the observations. 

Arguably, the field of nuclear astrophysics started with the monumental work of Hans Bethe, Subrahmanyan 
Chandrasekhar, and other great pillars of science. A humorous event might have happened (nobody really knows) in 
which during the 1930’s Hans Bethe told his girlfriend under a bright night sky watch that he was the only person on 
Earth who knew how stars shine. His girlfriend was apparently not impressed with his statement because girlfriends 
under a bright sky are more interested in something else. What Bethe knew and possibly nobody else did was that one 
gains a lot of energy by nuclear fusion of 4 protons into a ^He (a-particle) nucleus. Two protons turn into neutrons 
by positron emission (beta-decay process), releasing about 28 MeV (^^ 4.5 x 10“^^ J) of energy per event. At the 
time Bethe proposed a mechanism by means of which a series of reactions involving carbon and nitrogen nuclei served 
as catalyzers for the production of the a-particle (the CN cycle). The catalyzers are recycled, and the process is 
analogous to the work carried out by enzymes in biological processes. The orig inal model proposed by Bethe was 
later extended to include oxygen nuclei and is now known as the CNO cycle [^. Many decades passed by until we 
fully understood that Bethe’s model is responsible for energy generation in massive, hot stars, only. For example, in 
contrast to Bethe’s model our Sun catalyzes a-particles by means of proton-induced fusion reactions, also known as 
the pp-chain [30l |. 

The tunnel effect is responsible for the long time scale of stellar burning. Without tunneling there would be no 
stars as we know them. It is not a coincidence that all main physics related to the understanding of how stars shine 
occurred in parallel with advances in Nuclear Physics. The more we learn about Nuclear Physics, the more we learn 
about stars and the universe. Sorry, stellar modelers, you have to wait for us! The first application of the tunnel 
effect occurred in Nuclear Physics to understand the process of a-emission by heavy nuclei. For example, ^^*^Po emits 
a-particles that tunnel through a Coulomb barrier to get free outside the nucleus. This was explained by Gamow, 
Gurney and Condon who calculated nuclear half-lives by a-decay (or a-emission) using first principles of quantum 
mechanics available during late I920’s [3l|, [s^. The inverse process of capture of a-particles or any other charged 
particle (such as proton capture) is also explained as an inverse tunneling process through a repulsive Coulomb barrier. 

The immediate conclusion in this initial story-telling about nuclei and stars is that there would be no life in our 
Universe, or even no Universe as we know it (or the little we know about it!), without quantum mechanics ruling 
the physics of the nuclear world. The relevance of nuclear physics for stellar evolution is deeply rooted in processes 
generating energy and forming heavier elements. But there is a caveat: nuclear physics is arguably one of the most 
difficult research areas of all sciences. This is because the nucleon-nucleon interaction is not so well known and also 
because nucleons are composite objects ruled by the physics of quarks and gluons, effectively manifested through 
other particles which mediate the nuclear force (i.e., the mesons). Moreover, nucleons in nuclei are not so many to 
allow for simplifications traditionally available to understand the physics of bulk objects. In other words, nuclear 
structure and nuclear reactions are very difficult to handle using the underlying theory of quantum mechanics. In 
more mundane words, nuclear physics research is not for everyone. Only the strong and stubborn scientists survive 
in nuclear physics research. We hope that this statement will only motivate graduate students to choose research in 
nuclear physics. That is because difficult problems are more challenging and motivating for real scientists. 

It is extremely likely that life would not exist without the presence of the elements we have identified in nature. 
The light elements were produced in nuclear reactions during the Big Bang and heavier ones after the formation of the 
first stars. Carbon is a crucial element for life, although it is not ruled out that other elements such as silicon might be 
the basis of life in other strange planetary environments . A crucial step in understanding the origins of life as we 
know was the discovery of the triple-alpha capture process when three alpha-particles join to form a carbon nucleus. 
It was early concluded from many experiments that there would be no possibility to form carbon as abundant as we 
observe without resorting to the triple-alpha process, which is part of a path to the formation of heavier elements. 
If the triple-alpha capture process would be non-resonant then the reaction rate would also be too small in a stellar 
environment. It was necessary to postulate the existence of a resonance in the continuum of the beryllium nucleus 
(unbound) and also in the carbon nucleus close its alpha emission threshold, more precisely at 7.65 MeV [s^. The 
resonance in beryllium allows for a short time interval for the collision with another alpha-particle to occur, therefore 
the reason for coining it as the triple-alpha process. The resonance in carbon would be responsible for the large value 
of the capture cross section necessary to explain the existence of “abundant” carbon atoms in the Universe. 

This purely theoretical hypothesis, due to Sir Fred Hoyle was later confirmed by William Fowler and collabo¬ 
rators with an experiment at the Kellogg Radiation Laboratory at Caltech . The resonance is known as the Hoyle 
state or the nuclear state of life. Some think of it as a support to the anthropic principle: we exist, therefore this state 
[ssHstI in carbon is necessary, or has to be preconceived. Many scientist dislike the use of this principle in science and 
even argue that there is no principle at all, including the Nobel laureate Steven Wienberg |36l |. Such statements are 
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not even wrong, as used to say another Nobel laureate Wolfgang Pauli. Believing in it is very much like an ornament 
fish’s thinking: there must be a God, otherwise who would throw food in my aquarium? It would be embarrassing to 
find out that human brains are only as far able to think as a fish. Despite such spiritually inspired nonsense, due to 
its immense relevance, the triple-alpha reaction still leads to genuine theoretical and experimental interest |37 h39I| . 


B. Gluons, Quarks and Nucleons 


In 1949, again Hans Bethe [i^ introduced the concept of scattering length and effective range to describe nucleon- 
nucleon scattering at low energies. To simplify his description, in the absence of the Coulomb interaction the scattering 
phase shift expansion for a partial wave I is given by 

cot Si{E) = - hr/—-!-■•• , (1) 

ai Z 

where E = {hk)'^/2^, gi is the reduced mass of the system, ai is the scattering length and r/ the effective range. 

Bethe had no idea of the existence of an inner nucleon structure. Neither did Hideki Yukawa when he introduced 
a simple model to explain nuclear forces based on pion exchange [4l|. His model was later extended to include other 
mesons such as the rho and omega mesons. Now that we know that nucleons are made of quarks and gluons, we also 
know that the meson-exchange theory of nuclear forces is an effective theory by means of which the quarks and gluons 
degrees of freedom are hidden. Such a meson-exchange approach is very much like Bethe’s theory of Effective Range 
Expansion (ERE), Eq. ([IJ, which describes many features of the nucleon-nucleon forces in terms of the scattering 
length and the effective range. They hide the nuisances of the nuclear interaction at very low energies. 

Since we know that quarks and gluons flourish within nucleons, we also discovered that the QCD Lagrangian [4^ 

Cqcd = - \GlnGr . ( 2 ) 

is the basic theory describing the inner works of the nucleon-nucleon interaction. In Eq. 7 ^^ are Dirac matrices, 
ipnix) is the quark held, depending on space and time, indexed by to, n, • • •, and TOg are quark masses. Moreover, 


G% = d^A:-d,.Al + gf 


abc Ab 


A^ 


( 3 ) 


where 5 is a coupling constant, are the structure constants of SU(3) and Af{x) are non-Abelian gluon helds, also 
functions of space and time, in the adjoint representation of the SU(3) gauge group, indexed by a, 6 , • • •. The Gf^, 
are the QCD analogous to the electromagnetic held strength tensor, of quantum electrodynamics [i^j . 

There has been a strong effort to obtain nuclear properties from hrst principles by solving Eq. ©■ The hrst step, 
namely, reproducing hadron ground state and resonances has been rather successfully achieved, especially for light 
hadrons, by means of, e.g., Lattice QCD (LQCD) The LQCD is formulated to solve Eq. ([2]) on a grid or 

lattice of points in space and time [^. There is still a long list of hadronic problems to be solved before a fully 
consistent set of results is achieved, also for heavy hadron masses. It more difficult to describe scattering states from 
LQCD ( 4 ^ than hadronic ground state properties. It will certainly take decades or many human generations to extend 
the LQCD calculations to describe ground and excited state properties of many nucleon systems. First attempts in 
this direction are just at their infancy stage [s^. 

Using LQCD to obtain many bound and unbound nuclear properties also seems unnecessary and far fetched. It is 
much like killing a fly with a cannon. The world “effective” is in the root of a better solution for the problem. One 
does not need to use all the complex non-linear features of QCD for low energy nuclear physics because such features 
arise mostly at short ranges which are rather unaccessible in low energy processes. Thus, one needs an effective theory 
to handle nucleon-nucleon and multi-nucleon interactions in the same way as Bethe’s ERE theory was developed for 
scattering theory with wave mechanics. In 1979 Steven Weinberg proposed the use of Effective Field Theories 
(EFT) for low energy processes. In nuclear physics it means something like replacing the Lagrangian ([5]) by 


Cqcd Ceft = [idt + ^ ^ i^-Go{N^Nf + ^ [n^N | 


( 4 ) 


where N are nucleon helds for the isospin doublet, N = {p,n), which are functions of spacetime, and -1- 

2 -I- In this expression, p, is a constant to give the second term a correct dimension for any D. The 

constants C/ are called Low Energy Constants (LEG) and carry the hidden information on the short-range physics. 
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i.e., they carry the complex information about quarks and gluons inside the nucleons. The simpler Lagrangian (|4]) 
carries all required symmetries from the underlying QCD Lagrangian and can lead to more precise results as those 
dots on the right hand side are replaced by more terms in the expansion. Some say that what matters in EFT is how 
one treats those dots. 

In Quantum Field Theory (QFT), given a Lagrangian one can solve scattering problems by “reading” the Feynman 
rules out of it. Then one follows the usual procedure of QFT with regularization of integrals, renormalization due to 
infinite sums, and all those seemingly boring stuff. It is a completely different approach than solving the Schrddinger 
equation with a nucleon-nucleon potential. Thus, the solution of nucleon-nucleon problems using Eq. (jT]) can be done 
with increasing levels of approximation, namely, to leading order (LO), next-to-leading order (NLO), next-to-next-to- 
leading order (NNLO, or N^LO), etc. Thus, the level of precision of the calculation depends on how far does one go 
with the approximations (the dots). And they are (in principle) under control. In contrast, there is no predictive path 
to improve the Schrddinger potential method because nobody knows how to theoretically improve a phenomenological 
potential. EFT in fact gives a precise prescription of how to improve calculations [52l - [^ . Applications of EFT to solar 
nuclear reactions and to other stellar evolution in general have been reported in numerous works (see, e.g., (dll - l^ . 
For example, calculations of radiative capture reactions using EFT and lattice gauge theory have been reported (see 
e.g.. Refs. [6j, 1^). Calculations of the triple-alpha reaction and the state of life based on chiral-EFT have also been 
the recent focus in Nuclear Physics research 

Nuclear interactions are also strongly modified by the medium and solving the nuclear many-body Schrddingier 
equation is a task that has taken more than 80 years and is still being developed. The theoretical understanding 
of nuclear properties based on the QCD Lagrangian is a slowly growing computational effort, despite the advent of 
fast computers [t^. It might require another 80 years of theoretical efforts. But the use of QFT, EFT, QCD, and 
renormalization group techniques for low energy nuclear physics has been useful to unify all potential models used 
in nuclear structure calculations. For example, a particular interaction coined as Viow-k is shown to parameterize 
a high-order chiral effective field theory two-nucleon force. In the theory, a cutoff dependence can be used as a 
tool to assess the error in the truncation of nuclear forces to two-nucleon interactions and introduce low-momentum 
three-nucleon forces [nHzi- 

It is interesting that the nuclear physics community almost forgot about Bethe’s ERE which was quickly adopted 
in other fields. Nuclear physicists in the 1950-1990 were mostly interested either in nuclear spectroscopy in low energy 
nuclear physics, or in the study of a quark-gluon phase transition which might be discovered in relativistic heavy 
ion central collisions [s^. This phase transition might have occurred during the Big Bang and a colder version of 
it might also occur in the core of neutron stars. The research efforts and funding feeding the relativistic heavy ion 
community almost extinguished the scientific interest in low energy nuclear physics. Concepts such as ERE were 
being long forgotten in nuclear physics, although thriving in other areas. However, in the last three decades or so, 
the ERE method resurrected in the nuclear physics literature and spilled out to other fields. It was in fact another 
nuclear physicist, Herman Feshbach from MIT, that developed the concept of closed and open channels in reaction 
theory and who also introduced the concept Feshbach resonances. Both of these concepts became cornerstone tools to 
study atomic physics phenomena and in particular were well suited for the physics studied with ion and atom traps 
^|. ERE is now a standard tool used to understand many features in optical lattices [s^, Bose-Finstein condensates 
83, , or the unitary Fermi gases (ssl - fs^ . 


C. The Nuclear Chart, Varieties of Elements 

Assuming that we fully understand the nuances of the nucleon-nucleon interactions and how it emerges from QCD, 
we still would like to know how the nearly 250 stable nuclei emerge from a consistent theory. Most of these nuclei 
are composed of even number of protons and even number of neutrons, namely, even-even nuclei (see, e.g., [l^l. But 
even-odd and odd-even nuclei are also abundant. On the other hand, the number of odd-odd stable nuclei can be 
counted in our hands, nominally: they are rare and amount to only 5. Among these, 4 are very light, ^H, ®Li, ^°B, and 
^^N. The fifth stable odd-odd nucleus is very heavy and exists in an. excited, isomeric state, ^®°"*Ta. It has a high 
spin, 9“, which makes its decay to the ground state, I"*", by y-emission or through /3-decay extremely unlikely since 
these decay processes only change the angular momentum by one unit. Thus to decay to the ground state, ^®°Ta, 
with spin and parity 1+, and whose half-life is about 8.15 h, there must be an exceedingly weak and thus very slow 
multi-photon emission or a similarly very slow /3-decay. As such, this isotope of Tantalum is stable in this isomeric 
state {E* = n.lkeV), and it is the only stable nucleus to exist in such a state. Another feature of ^^^'"Ta is that it is 
the rarest isotope in nature amounting to only 0.012% abundance in the natural ore composed mostly of the odd-even 
nucleus ^^^Ta (99.988%) (H. 

Isotopes of neutron-rich nuclei extend till the neutron drip line, where the neutron separation energy is zero, while 
isotopes of proton-rich nuclei are fewer and extend to the proton drip-line which lies closer to the valley of the 
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stability than the neutron drip line. The drip lines act as natural boundaries of bound nuclei. One important fact to 
mention is the nonexistence of stable nuclei with A = 5 and 8. This fact has an important consequence on Big Bang 
nucleosynthesis and especially on the carbon production. Recently it has been possible to produce neutron-rich and 
proton-rich radioactive nuclei as secondary beams, which allowed researchers to extend the study of nuclear structure 
to these exotic species of nuclei BlMl. This area of research has become a priority in many countries and major 
projects involving the construction of large facilities of accelerators are in the process of being developed in Germany, 
the USA, China, Japan and elsewhere. In Brazil, the activity in this area is centered in RIBRAS, installed in the 
Nuclear Physics Laboratory of the Institute of Physics of the University of Sao Paulo [H, 

Another important extension of nuclear research has been in the production of the so-called superheavy elements, 
with proton numbers exceeding 100. These elements are produced through the fusion of a medium mass nucleus with 
a heavy one resulting in an excited compound nucleus that, after emitting several y-rays and cooling down, starts 
emitting a-particles till reaching the superheavy element that survives enough time to be studied both physically and 
chemically (see, e.g.. [o^). Some of these elements are: Z = 101, Mendelivium Md, Z = 102, Nobelium, No, Z = 103 
Lawrencium, Lr, Z = 104, Rutherfordium, Rf, Z = 105, Dubnium Db, Z = 1 06, S eaborgium, Sg, Z = 107, Bohrium, 


Bh, Z = 108, Hassium, Hs, 109, Meitnerium, Mt, Z = 110, Z = 112, etc. 95l - ll03j| . 


There is a huge effort in Nuclear Physics to understand how nuclear masses emerge from a consistent theory of the 
nuclear forces. For medium and heavy mass nuclei, the mean field theories are the tool of choice. Starting from a 
nucleon-nucleon interaction including two and three-body terms, such as 


= E 

i<j 


,( 2 ) 


E 

i<j<k 


(3) 

%k 


( 5 ) 


the energy or nuclear masses can be calculated by means of m 


^=E 


2m N 


I E 1^12 I ^ E kl23| , 


( 6 ) 




i,j,k 


where m^r is the nucleon mass, and A mean s an ti-symmetrization of the matrix element. If a Skyrme contact-like 
interaction is used for v, it has been shown 104| that the above functional can be easily calculated if the nucleon 
wavefunctions |i) are obtained from an iterative procedure such as the Hartree-Fock method. Medium corrections of 
the nucleon-nucleon interaction by means of the Brueckner method or the inclusion of pairing by means of a BCS 
or the Bogoliubov metho d, ca.n im prove calculations of the functional 0. Inclusion of relativistic ef fects have als o 
shown to be of relevance [ll2lll22j |. Several groups are pursuing these calculations, e.g., those in Refs, 
and many open questions remain as, e.g. how well one can treat parin g in nucl ear matter. The Beijing group has 
made many relevant contributions to the progress in this research topic |ll3l4l2ll |. 

One of the most anxiously expected developments in this field is the use of time-dependent many-body techniques 
to unravel real-time excitation, decay, transfer and fragmentation of nuclei during nuclear collisions. It is a very hard 
task which requires the use of supercomputers. Advanced work on this subject has been reported in Refs. [l24l - ll28 |. 


D. Nuclear Sizes 

Owing to the great amount of energy required to excite or change nuclear structure, as compared to atoms, it 
is difficult to use external fields such as electric or magnetic fields to change the properties of nuclei or the basic 
nucleon-nucleon interaction. This is to be contrasted to the case of atoms, where the atom-atom interaction can be 
altered through the application of an appropriate magnetic field which affects the so-called Feshbach resonances. With 
this effect the atom-atom interaction can be changed from repulsive to attractive, and going through zero, making 
possible the study of cold gases under these different conditions. In nuclei no such liberty is available and the only 
way available for inflicting changes is the direct interaction with other nucleons (capture reactions) or other type 
of hadrons. By adding more neutrons or protons, new types of nuclei are produced as the case in nucleosynthesis. 
Some of these changes have become possible with the advent of the field of nuclear fragmentation in which nuclei 
are broken up and fragments rich in neutrons or protons are separated and further accelerated as secondary beams, 
whose properties are then discerned and analyzed. Nuclei, such as ^^Li, have been shown to have an rms radius as 
large as that of the heavy nucleus of lead, ^°®Pb! Thus, larger nuclei with special prope rties such as the existence of 
a halo of excess neutrons or protons (such as ^^Be, ^^Li, and ®B) can now be produced 

One of the important features of nuclear structure is the fact that nucleons are identical fermions and are thus 
subject to the Pauli exclusion principle. Excess nucleons can only occupy the outer most orbits. Such restriction is 
removed if the added particles are other types of hadrons. One such example is the hyperon A which is a strange 
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nucleon whith a strangeness quantum number of 1. Once the A is introduc ed in to the nucleus, it does not suffer the 
effect of the exclusion principle and can migrate to the center of the nucleus 13l| . Through the attraction it exerts on 
the other nucleons, these tend to aggregate closer to the center, making the normal nucleus shrink in size. In the hyper 
nucleus ^Li ( = A + ®Li), the radius of ®Li is 20% smaller than that in free space. It would be interesting to study 
such shrinkage in siz e in heavy nuclei, such as ^°®Pb, through the production and investigation of the superheavy 


hyper-nucleus f’^Pb [13 

Other ways to influence the nuclear structure have been suggested. Very recently, multi-MeV zepto-second coherent 
laser pulse by backward Compton scattering of optical laser light on a sheet of relativistic electrons has been proposed 
to provide a huge number of MeV photons into the nucleus w hich could either creat e a plasma of nucleons or to 
excite collective modes. These lasers are being developed at ELI [I33j| and IZEST 1341. These developments in laser 
technology will enable nuclear researchers to go beyond the usual studies of the nuclear response, and to actually 
inflict changes in the properties of nuclei which are hitherto not possible and could reach the realm of freedom of 
manipulation practiced in cold atomic gases and Bose-Einstein condensation research. 


III. NUCLEAR REACTIONS 
A. The Optical Model and Beyond 

Research in nuclear reactions accompanied closely the development of models of nuclear structure [l35j| . Two aspects 
of nuclear structure were in apparent conflict. The shell model which assumes a mean field felt by the nucleons and 
a corresponding long mean free path, and the Bohr model of the compound nucleus which asserts that the nucleons, 
once captured suffer many collisions with the other nucleons and with the mean field, which implies a short mean free 
path. These conflicting models were later shown to be different manifestation of the same underlying nucleon-nucleon 
force and the Pauli exclusion principle. A model for the nuclear reaction which is an extension of the shell model 
to positive energies, was then developed, where the nucleus is taken to exhibit refraction due to the real mean field 
and diffraction due to the absorption caused by the formation of the compound nucleus, was developed by Eeshbach, 
Porter and Weisskopf |l36j | , and is based on the use of a complex average potential felt by the impinging nucleon. 
The optical model equation is 

(A:-ht7)|«'(+)) =E|^'(+)) . (7) 

where K is the kinetic energy operator and U{v) = V{r) — iW{r) is the complex potential, whose real part, V{r) is 
related to the mean field used in the shell model (at negative energies), while the imaginary part, — IT(r), accounts 
for the flux lost into the formation of the compound nucleus which is treated separately using the statistical model. 
In fact, a simple manipulation of the above Schrodingier equation can be used to derive a continuity equation through 
which an absorption cross section can be derived, 

<7absiE) = ^(4-(+)|tT(r)|^'W) . (8) 

The above cross section accounts for the compound nucleus, treated within the optical model, as a sink of flux. 
This cross section can be expanded into partial waves, yielding, 

OO 

<JaUE) = ^Y.^2l + l)Ti{E) , (9) 

1=0 

where the Ith transmission coefficient, Ti{E) is 

A h pco 

Ti{E) = -j^ dr\i;l+\r)\^W{r) , (10) 

and is the radial wave function of the optical model with scattering boundary conditions. 

In the case of fusion cross sections, Eq (l9|) can be used with X = 1/k = /2mE = h/mv being the reduced 

wavelength (remember the reduced Compton wave length Xc = h/mc). Eq. ([9|) can be interpreted as if the cross 
section is proportional to ttA^, the area of the quantum wave. Classically, different parts of the wave have different 
impact parameters and different fusion probabilities, T;. Large impact parameters correspond to large angular mo¬ 
menta, leading to the weight {21 -|- 1). Eor fusion, it is very useful to use the concept of astrophysical S-factor, such 
that 

(^f{E) = ^S{E) exp [-27r?7(£')] , 


(11) 





where r]{E) = ZiZ 2 e^/hv, with v equal to the relative velocity. The exponential approximately accounts for the 
barr ier t ransmission. While cr{E) decreases rapidly as the energy decreases, the astrophysical S-factor remains rather 
flat [T^ . 

Empirical optical potentials used in the analysis of elastic scattering of nucleons and heavy ions have been developed. 
Pion-nucleus scattering was also extensively studied using semi-microscopic optical potential based on salient features 
of the pion-nucleon scattering t-matrix including the excitation of the A-resonance. In the case of heavy-ions, we 
mention the recent development of the Sao Paulo Potential, based on the double folding of the dens ities of th e two 
nuclei in conjunction with an effective N-N interaction including non-locality owing to Pauli exchange [niiiii. This 
potential has been quite successful in accounting for the elastic scattering of both stable and radioactive nuclei, after 
the a ddition appropriate polarization potentials t hat simulate the effect of the couplings to strongly coupled channels 
[l40j| . such as the breakup one in the latter case [l4l| . Its use in coupled channels calculations has bee n pro ven to be 
quite successful in dealing with the scattering of nuclei from deformed targets, and in fusion reactions [l42|. 


B. Compound, Pre-equilibrium and Direct Reactions 

Low energy reactions are dominated by the mechanism of the formation and subsequent decay of the Compound 
Nucleus (CN), which is the nucleus formed from the capture of the projectile by the target nucleus. This compound 
system is formed at a relatively high excitation energy and angular momentum. According to Bohr’s hypothesis, the 
formation and decay of the CN are independent as the whole process is statistical and it takes time for the system 
to equilibrate after the capture process. Accordingly it is assumed that the cross section to go from channel a (the 
entrance channel) to a final decay channel P, is for a given value of the angular momentum, 

CTa./3 = -^(2^ -b l)?7a?7/3 • (12) 

Summing over P gives the total absorption cross section for angular momentum I, ^a,p = (tt/ k'^)r]a, V/s which 
can be further simplified by recognizing that the transmission from channel a, = iP^/ tt{21 + l))crct, can be written 
as T„ = ria Vp- Thus Tq, = (^^ ■ These manipulations leads to rja = Ta/ ^( 3 , and the cross section 

becomes, 

ac.,p = ^i2lEl)^^^. (13) 

The above expression is the Hauser-Feshbach ( HF) cross section [l43j| . The importance of the work of |l43l| which 
is an extension of the Ewing-Weisskopf theory |144| . is the observation that the transmission coefficients Ts, where 
S designates any decay channel of the compou nd n ucleus, including the entrance channel, is directly related to the 
absorption cross section of the optical model |144| . The case of compound elastic, a = /3, is distinct as there are 
obvious corre latio ns between the entrance and exit channels, and the HF cross section above has to be multiplied by 
a factor of 2 |l45j| . The angular distribution of compound nucleus cross section is symmetric around 0 = tt, and in 
fact is isotropic as expected of a statistical process. 

Higher energy data on particle spectra have indicated deviations from the pure compound nucleus evaporation form, 
indicating the operation of another mechanism of decay. This particle emission process is associated with the excitation 
of 2p-lh (2 particle - 1 hole), 3p-2h, etc, configurations in the compound nucleus before a fully equilibrated system is 
reached. Such pre-equilibrium emission has been the subject of investigation both experimentally and theoretically. 
The angular distribution of the emitted particles is, however, not isotropic, indicating the contribution of yet another 
mechanism besides that of the CN. This other mechanism involves the excitation of the target nucleus through the 
excitation of Ip-lh, 2p-2h configurations with the projectile particle remaining in the continuum. The combined effect 
of multi step compound and multistep direct processes constitutes the model of pre-equilibri um r eactions which is used 
in Ref. [l46j| . For a recent review of compound nuclear pre-equilibrium reactions see Ref. [l47j| . 


C. Feshbach’s Theory 

Nuclear reaction data indicate clearly the existence of two types of processes. The ones dominated by the compound 
nucleus which are slow processes and treatable with the statistical model of Bohr and formalized through the Hauser- 
Feshbach theory, and others, fast processes, describable by the optical model equation as generalized to many channels. 
Processes, such as inelastic scattering, transfer reactions, and breakup reactions are commonly called direct reactions 
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and characterized by forward peaked angular distributions with oscillations indicative of their coherent, non-statistical 
nature. To accommodate both types of reactions, the slow, compound ones and the fast, direct ones, Feshbach 
[148L Il49l |. developed a formal theory based on projection operators. 

Call the open channels space projection operator, P, and the closed channels space Q, with P+Q = 1 and PQ = QP 
= 0, PP = P, QQ = Q, then the Schrddingier equation of the colliding nuclear system can be written as, 

{E - PHP)P\'l’) = PHQ\'l’) (14) 

{E-QHQ)Q\'i')=QHP\'i/) . (15) 

These equations summarize the whole theory of nuclear reactions. To describe the direct reactions one eliminates the 
compound nucleus Q-space and averages out the corresponding resonances, to obtain the following effective equation, 

= ( 16 ) 

The Hamiltonian H in the above equations is the sum H = K + hi + h 2 + V 1 . 2 , where K is the kinetic energy 
operator, hi is the intrinsic hamiltonian of nucleus i, and Fi ,2 is the real interaction operator between the two nuclei. 
The effective polarization operator 


PHQ 


{E - QHQ) 


QHP = PHQ 


{E - QHQ + il) 


QHP , 


where the energy averaging interval, I, is much larger than the average width of a resonance. This equation is complex 
and by adding it to PVi^ 2 P, defines the complex optical potential operator. 


Veff =PVi,2P +PHQ 


- QHP, 

E-QHQ + il^ 


(17) 


and the optical model equation becomes. 


{E-K-hi-h2-Veff)P\'i'> = 0 . (18) 


This equation is an equation for the optical wave function > which contains many components (channels). As 
such, it is a set of coupled channels equations that describe the direct reactions involving the channels projected by 
P. The compound nucleus is now completely hidden in the complexity of 14//- To extract the contribution of the 
compound nucleus, one needs to resort to statistical considerations involving the fluctuation component of >. 
The treatment of this issue relies on the use of assumed random properties of the P — Q couplings and the whole 
apparatus of quantu m chaoti c scattering theory is employed. The Hauser-Feshbach cross section and the Ericson 
correlation function |l5fll Il5l| are then obtained, making the above theory a complete one. 

The set of equations Eq. (US, constitutes the Coupled Channels T heory (CCT) of nuclear reactions. Most reactions 
at low and intermediate energies require treatment with the CCT [l52j| . At higher energies or weak couplings, a 
perturbative treatment is adequate, through the use of the Distorted Wave Born Approximation (DWBA). This 
theory has been very useful in the study of the nuclear structure as the amplitude is linear in the coupling that 
induces the transition and accordingly spectroscopic information are unambiguously extracted Q. Variance of this 
theory, which is commonly used in the case of transfer of nucleons in cases involving deformed nuclei, is the Coupled 
Channel Born Approximation (CCBA), where the distorted waves are generated from a coupled channels calculation, 
whereas the transfer coupling is treated perturbatively. Such theo ry is still in use in the case of the scattering of 
exotic neutron- and proton-rich nuclei from deformed target nuclei |15,8l |. 


D. Coupled Channels 

The advent of radioactive beams has brought into focus several important features of reaction dynamics. The CCT 
is an example. The wave function P|4'(+) > is taken to be composed of a sum of several terms representing the 
important channels operating in the collision of a stable, or neutron- or proton-rich projectile with a target nucleus. 
Writing P = Pf,i + -|- Ptrans + P' renders the optical model equation a set of coupled equations for the elastic 

channel, Pg/, the breakup channel, P{,„p, the transfer channel, Ptrans, and the other channels that are treated in 
an average way, P'. For exotic halo nuclei, such as the one-neutron halo ^^Be, the two-neutron halo ^^Li, ^°C, or 
one-proton halo, ®B, the breakup coupling is quite important owing to the close to threshold breakup channel. As 
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such, it becomes important to consider Pei and Pbup in the projected P-space optical equation, when treating elastic 
scattering. 

The adiabatic model, i.e., the neglect of the breakup Q-value, is generally used at higher energies, as done in 
and [l55j| . In |l54j the Glauber approximation while in |l55|| the usual adiabatic Sc hrdin ger model are used for the 
radial wave function. The result of the analysis of the total reaction cross section in |l54| of ^^Li, ^^Be, and ®B and 
a target, demonstrated that the breakup channel coupling necessarily increases the reaction cross section and 
accordingly requires a larger rms radius of the halo nucleus, as seemed to be required by the data on the interaction 
cross section (the total reaction cross section without the proj ectile inela stic or breakup contributions) at the time, 
which were analyzed using the optical, eikonal, approximation This result is easily understood using the 

general structure of the transmission coefficient Ti(x) = T{b,x) = 1 — exp [—2y(6, x)], where x(&, x) is the imaginary 
part of the eikonal phase, b refers to the impact parameter and the parameter x refers to the se parat ion distance 
between the core and the center of mass of the excess nucleons. The original, Tanihata analysis |l57j |. considered 
the optica l average exp [—2 (x(&, x)),,,] whereas the adiabatic model deals with (exp [—2x(&, x)])^. Jensen’s inequality 
[l58l Il59l | comes into play to dictate that (exp[— 2 x]) 3 , > exp [—2 (x)^,], and thus the transmission coefficient in the 
adiabatic model Tad{b) = 1 —(exp [—2^)^ is smaller than the optical transmission coefficient, Topt{b) = 1—exp [—2 (x)^;] 
rendering the reaction cross section aa^ad = Stt bdbTad(b), smaller than the optical one, (jR^opt = 27r bdbTopt(b). 
To compensate for the reduction in (TR_ad compared to (iR^opt, a larger rms radius, e.g. for ^^Li, < rh > = 3.55 fm, 
than the one extracted by Tanihata 
the treatment of the elastic scatterinj 


|157| . 3.1 fm, is requ ired to fit the data.The same considerations were applied in 
of ^^Be + Il55|. Im provem ent of the adiabatic mo del was t hen developed by 


the Pisa [I60l - ll67j |. Commerce |168L Il69l |. Brussels |l70l - ll72j| , Osaka and Seville groups, through 

the relativistic Continuum Discretized Coupled Channels (CDCC) model and the dynamical eikonal model. 

At low energies the fusion process becomes important. A lot of attention has been dedicated to the influence of the 
excess nucleons on the tunneling probability which dictates the value of the fusion cross section at near Coulomb barrier 
energies. The discussion of this process and other reaction processes at low energies requires the solution of coupled 
channels equations involving the projections Pei, Pbup, and Ptransf exactly. The breakup channel involves three or 
four clusters in the continuum, requiring for its treatment a formidable three- or four-body scattering calculation. 
In practice, however, the contin uum is di scretized into a finite number of pseudo-states treated as inelastic channels. 
The resulting CDCC equations, [UMli, are solved for a truncated number of pseudo states using known numerical 


methods. Recent improvement and extensions of the CDCC have been made |184 Il85| . Computationally the CDCC 
relies on an ad hoc truncation of the discretized continuum by keeping only a few pseudo-states. The effect of the 
neglected pseudo-states on the convergence of the results is se ldom analyzed. Recently, the idea of using statistical 
methods to treat these neglected pseudo-states was advanced |l86l| . and a model for a Stat istieal CDCC (sCDCC), 
based on Quantum Chaotic Scattering Theory (see below) is currently being developed [l87j . 

The results of CDCC calculations for fusion reactions have taught us several things. The excess nucleons in the 
halo nuclei results in two important effects. The excit ation of low lying dipole mode, the Pigmy Dipole Resonanee 
(PDR) (also known as the Ikeda resonance) [HMn, and the close-to-threshold breakup. These two effects were 
found to influence the fusion in two distinct ways. The breakup coupling results in a suppression of fusion at energies 
above the top of th e Coulomb barrier, whereas the dipole excitation results in an enhancement of fusion at sub-barrier 


energies |202l l203j . This latter effect is intimately related to the longer extension of the matter density of the halo 
nucleus. An extensive effort, both experimental and theo retical, has been dedicated to the low energy reactions of 
exotic nuclei and several reviews were published, [2041 - 1207 1. 

The interest in fusion of neutron-rich nuclei stems from the desire to produce a compound system that could survive 
long enough to be studied and analyzed. Heavier nuclei with large amount of excess neutrons have been produced 
as secondary beams, and the existence of the Pigmy resonance was clearly established in, e.g., ^^"‘Sn |l89l . Il9ll - ll93 |. 
These nuclei with a thick neutron skin, were produced at high energies. It would be very interesting to produce the 
neutron skin nuclei at lower energies and study their fusion with a heavy target to test the idea of an enhanced fusion 
and a potential production of a superheavy nucleus with a Z and A beyond the known ones. Such investigation will 
become possible in the near future with the new nuclear facilities being built around the world. 

It is worthwhile mentioning that recent experiments have also investigated the nuclear neutron skin by means of 
completely different physics me thod s . Fo r example, we cite the ^°®Pb Radius Experiment (PREX) based on parity 
violation in electron scattering Such ex perim ents seernto be l ess dependent o n the data interpretation 

base on a model description for scattering by pions [l96| | , protons [I97l4l99| , or antiprotons [loSlloIl. 


E. Incomplete Fusion 


Besides complete fusion |208l |. which involves the capture of the whole projectile, the breakup may lead to the 
capture of one of the fragments, resulting in what is known as incomplete fusion. Other names of this process are 
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used in the literature, such as massive transfer, inclusive non-elastic breakup, surrogate reaction, etc. This process is 
important at it supplies a mean to study the fusion of a nucleus which is otherwise difficult to produce as a beam, 
such as neutrons (in a deuteron induced reaction). Several methods and theor ies have been proposed to calculate the 
incomplete fusion reactions. We mention the ones which are currently in use |209l - l214 |. 

The standard treatment rely on the spectator model, which says that the two-cluster projectile a = x + b interacts 
with the target only through the participant nucleus x, leaving the spectator b unaffected and only suffer elastic 
scattering. The cross section for observing b is invariably derived to be, 


(Pa 

dfltdEb 


hVa 


> 


(19) 


where Eb is the energy of the outgoing fragment b, and pb{Eb) = rubkb/{8 tt^ h^), is the density of states of b. The 
source function, pxijx), is the wave function of the participant fragment x in the projectile as it reaches the target 
nucleus. T he irnaginary p art of the x-A optical potential is designated by 
In Refs. |21 CiL |211L l214l |. the source function is the overlap 


/3x(rx) = ^(rb) xi+)(r„,rb)(/)a(rb-r^)), 


( 20 ) 


where x is the optical model wave function (distorted wave), and (j)a is the intrinsi c wave fun ction of the projectile. 


The coupling interaction in Refs. [209142llj was taken to be I4_b (post), while in [212l 


the inte raction is the 
2091 . l212l . |213|| , the incident 


difference in the optical potentials, Ux{xx) + Ub{rb) — Ua{xa) (prior). Further, in Refs. 

wave function of a is taken to be beyond the distorted wave. A debate has been going on conc erning wtiicti ot tJie aPove 
approaches is more appropriate for the calculation of incomplete fusion |210l . l21ll l215l l216j| . This debate continues. 

The im portant po int to mention here is that the structure of the cross section above is similar to the one suggested 
by Serber [2171 l218j , namely the spectrum of the spectator fragment is proportional to the squared Fourier transform 
of the intrinsic wave function of the proj ectile , tim es the total reaction cross section of the participant fragment. 
Further analysis was performed in Refs. [2id . l21l| using the eikonal approximation for the distorted waves gave 
for the incomplete fusion cross section the simple form of an integral over impact parameter of terms of the type, 
((1 — Tb)Tx), where T is the optical transmission coefficient, and the average is over the internal motio n of b inside 
the projectile. If the Tb is set to zero and the average is ignored one recovers the Serber expression |2I 


d?a 

dflbdEb 


Pb{Eb)\(l>a{^b)\'^crx{Ex) ■ 


( 21 ) 


The precise determination of the incomplete fusion is important not only for the purpose of capture reactions of 
only a part of the projectile, such as neutrons in deuteron induced reactions, but also in fusion studies as in many 
instances the data give the total fusion which is the sum of the complete fusion plus the incomplete fusion. Therefore 
to get the complete fusion one needs to subtract from the data a believable incomplete fusion cross section. The 
research in the area of incomplete fusion is currently pursued by several groups. A semiclassical model has recently 
been developed 12191 1. In fact only in this year three papers have been written and will appear soon published in the 
journals [2^222l| . 


F. Quantum Chaotic Scattering 


A subject of continues inter est in nuc lear physics research is the understanding of fluctuations in the cross sections. 
Back in the early 50’s Wigner [223l |224| recognized that it is meaningless to try to analyze the many resonances seen 
in the neutron capture cross sections in the region of overlapping resonances. He suggested that a statistical treatment 
is more appropriate and introduced the concept of random matrices for the purpose. The idea is to assume that the 
nuclear Hamiltonian be considered a member of an ensemble of random matrices and averages of different quantities 
be perfumed through ensemble averages. Th e Ra ndom Matrix Theory (RMT) of nuclear reactions was born. The 
RMT was then greatly developed by Dyson [225[| and application to fluctuation phenomena in other fields besides 
nuclear physics ensued. 

At the same time Ericson [Miiiliig introduced the cross section correlation function as a measure of the 
degree of coherence in the othe rwise chaotic nuclear system. Through an analysis of the correlation function or the 
average density of maxima [230[| one is able to extract the correlation width, which had been previously estimated by 
Weisskopf to be 


F 


corr 


D 

W) 



p 


( 22 ) 
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where D is the average spacing between the overlapping resonances. These developments were quite important as 
they are universal in nature and can be employed in many systems exhibiting fluctuations in the observables. This 
syner gy, s upplied by research in nuclear physics, attests to the richness of the ideas and concepts introduced in the 
field ( 23 1| . As a matter fact, recent applicati on o f RMT has been mostly in mesoscopic systems, such as electronic 
conductance in quant um dots a nd graphene |232| . The test of RMT has been made possible using, among others. 


microwaye resonators [233l l234| . 

The Quantum Chaotic Scattering Theory relies basically on an expression of the S'-matrix which exhibits its relation 
to the random Hamiltonian of the system. It is normally written as, 


5'(£:) = 1 - 2TTiW\£ -H + iirWW'^y'^W 


(23) 


In the equation aboye H is the random Hamiltonian which pertains to one of the universal classes of random 
matrix ensembles, the Gaussian Orthogonal Ensemble^ GOE, the Gaussian Unitary Ensemble, GUE or the Ga ussian 
Simpletic Ensemble, GSU, considered by Dyson and excellently reviewed by Bohigas and Giannoni |225l . l226l | . The 
coupling matrix W couples the internal degrees of freedom to the open channels, and is taken to be fixed (not random). 
The observables such as the average cross section and the correlation function are then calculated by performing an 
average over the ense mble to w hich iJ pertains of products of 5'-matrices of the type given above. For more details 
we refer the reader to |233l l234j . 

Experimentally, the emergence of chaos in quantum systems can be verified by looking at the statistical properties 
and eigenvalues of the wave functions, and in in the fluctuation properties of the scattering matrix elements of some 
systems, such as quantum dots, quantum wires and Dirac quantum dots (quantum dots on a graphene flake, where 
the electrons are massless and obey the Dirac equation rather than the Schrdinger equation), etc. Microwave billiards 
provide very useful systems to study quantum chaos, because they already contain a degree of chaoticity in their 
classical dynamics. In f act, the e igenvalues and wave functions of quantum microwave billiards have been studied by 
Achim Richter’s group [227l - l229l| with a high precision insight into quantum chaos phenomena. 


IV. NEW NUCLEI 


A. Halo Nuclei, Efimov States, and Borromean Nuclei 


A rather simple but ingenious experiment reported in 1985 by Isao Tanihata and his grou p on interaction cross 
sections of light nuclei close to the drip line was the seed of a new era in nuclear physics |I56l| . This experiment and 
others following it, have shown that some nuclei such as ^^Li possess a long tail neutron distribution. The long tail 
is due to the low binding energy of the valence nucleons. It is simple feature, but one that nobody saw before, and 
it took Tanihata’s genius to re alize the apparently obvious nuclear property. His experiments were put within a nice 
context by Hansen and Jonson [235j| who relied on t he analysis of Goulomb breakup of ^^Li, shown to be enormous due 
to the loosely bound character of the nucleus [236j| . After Hansen and Jonson’s paper was published the community 
suddenly realized Tanihata’s discovery and the number of works and citations to the emerging field of exotic, halo, 
and neutron-rich nuclei skyrocketed. 

These developments lead to the then existing nuclear facilities in Europe (GSI/Germany and GANIL/France), in 
the USA (NSCL/MSU), and in Japan (RIKEN) to almost entirely turn their beams and detectors to the production 
of radioactive beams. Experiments with nuclei far from the stability became a routine in nuclear physics. The opening 
of this new research field in nuclear physics was one of the reasons why the U.S. National Science Foundation decided 
not to close the NSGL facility in Michigan State University (MSU) in 1992. The Indiana Gyclotron facility was not so 
lucky at the time, perhaps because they missed those historical developments. Since then new facilities have evolved 
and planned in all continents, in particular we mention the n ew RIKE N radioactive beam facility and the ERIE/MSU 
and FAIR/Darmstadt facilities presently under construction 237Ll238 |. In particular, it took the relentless and diligent 
work of Konrad Gelbke and collaborators to elaborate and get funds for the FRIB/MSU facility in the U.S. Also in 
China and Korea other facilities are being planned or under construction. For a r eview of the theoretical efforts 
conducted at the period of 1990’s, we suggest the reader to read Refs. gsmiii. 

With so much of science at stake it was imperative to obtain funding for new nuclear physics laboratories with 
the sole goal of studying nuclei far from the stability line. It is worthwhile mentioning that before this era, most 
physics studied in those facilities had to do with central coll ision s with the purpose to study the equation of state 
(EOS) of nuclear matter at high densities and temperatures [242j| . This EOS (pressure versus density), specially at 
low temperatures, is crucial for understanding the physics of supernova explosions and neutron stars. In fact, masses 
and radii of neutron stars are constrained by the EOS of nuclear matter. In contrast to the study of EOS, the physics 
of radioactive secondary beams was mostly driven by peripheral, direct reactions. Stripping, Coulomb and nuclear 
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excitation, two- and three-body breakup, and other direct reactions were and remain the principal tools to access 
spectroscopic information of interest to model and develop theories for short-lived nuclei, far from the stability line. 
Their relevance for astrophysics purposes is discussed below. 

Typically, neutron binding energies are of the order of 8 MeV for most stable nuclei. But neutron separation 
energies as low as few keV have been observed in short-lived nuclei. The long (visible only in a logarithmic scale!) tail 
of the weakly-bound neutrons and their corresponding densities is known as “halo” and the nuclei with this property 
are known as halo nuclei. Not only the small separation energies, but also nucleon-nucleon correlations are vital for 
the halo properties. As a typical example, neither ^°Li nor two neutrons form a bound system, but together in ^^Li 
they bind by about 300 keV [s^. The nucleus ^^Li became a sort of nuclear superstar in the 1990’s. There was not a 
single month or weeks without a nuclear physics preprint reporting either an experiment or a theoretical calculation 
on ^^Li. Nuclear physicists also love fashion, fame, and applauses, not differently than anybody else. 

Three-body systems theorists were very happy to see that ^^Li fell in their category of delicate structur es ar ising 
from loosely bound three-body systems. In fact, ^^Li is a prototype of a structure known as the Efimov effect 
the appearance of (many) bound states in a system of two-body subsystems with large scattering lengths [24C 
Jan Vaagen named such systems as Borromean system s, be cause they reminded him of the heraldic symbol of the 
aristocratic Borromeo family from the fifteen century [239|. In fact, the Borromean rings shown in the heraldic 
symbols are three inter-connected rings, such that if one is cut loose, the remaining two also become free. A large 
activity in the area of three and few-body physics have profited from this early work on nucle y ph ysics in the 1 980’s 
and 1990’s. Efimov states in Borromean systems have become a common feature in atomic [246j | and nuclear 2471 
physics. 


B. New Magic Numbers, Clusters, Majorana Particles 


In 1949 Maria Goeppert-Mayer and Hans Jensen proposed that nucleon forces should include a spin-orbit interaction 
able to reproduce the nuclear magic numbers 2, 8, 20, 28, 50, 82, 126, clearly visible from the systematic study of 
nuclear masses fli- Now we know that as nuclei move away from the line of stability {Z ^ N), where Z (N) 
is the proton (neutron) number) those m agic numb ers might change due to correlations and details of the nucleon- 
nucleon interaction, e.g., the tensor-force |248l - l254 |. These and other properties of nuclei far from the stability and 
in particular those close to the drip-line (i.e., when adding one more nucleon leads to an unbound nucleus) increased 
enormously the inter est for low-energy nuclear physics, in particular for the prospects of its application to other 
areas of science [255j| . The applications in astrophysics are evident. For example, the rapid neutron capture process 
(rp-process) involves nuclei far from the stability valley, many of which are poorly, or completely, unknown. 

It is not easy to predict what one can do with short-lived nuclei on Earth. Sometimes, they have lifetimes of 
milliseconds, or less. It is tricky to explain to the laymen the importance of such discoveries in the recent history of 
nuclear physics. For example, in 1994, the unstable doubl y-ma gic nucleus ^°°Sn was discovered at the Gesellschaft 
fiir Schwerionenforschung (GSI), Darmstadt, in Germany [256j| . At the same time the laboratory developed a new 
cancer therapy facility based on the stopping of high energy protons in cancer cells. To celebrate the new facility, 
the laboratory invited high-ranked authorities from the German government. They surely brought along the main 
news media in the country. Wisely, the laboratory management included the discovery of ^°°Sn alongside the proton 
therapy facility in the celebration agenda. At one moment a reporter asked a nuclear physicist what is ^°°Sn good 
for practical purposes, what was followed by a long silence and a curious answer: “I think one can develop better, 
lighter, cans.” Maybe so, but such a can would last fractions of second! The only sure prediction for applications of 
short-lived nuclei is that it is impossible to understand stellar evolution without a dedicated study of their nuclear 
properties. 

Clustering phenomena in nuclei is a difficult problem, which requires physics beyond the shell model (SM). The 
SM describes protons and neutrons moving in shells. Other models start from a molecular viewpoint that employs 
composite particles, such as the a-particles as the building blocks for nuclear structure. A major task is underway to 
connect the two seem ingly disjoint points of view and how nuclear excitations can emerge from a unified picture of 
the nucleus 2571 - 1272 1. 

Majorana had already shown in 1937 that if the neutrino would be its own anti-particle then the beta-decay theory 
would not change. If the neutrinos are their own anti-particle (Majorana neutrinos), then neutrinoless beta-decay 
in double beta-decay is also possible, with two neutrons being converted into two protons in the decay. But this 
is a very rare process and extremely suppressed, what makes it difficult to access experimentally. There has been 
some early hints that the proc ess might occur in nature but up to the present date there has been no definitive 
proof of its existence |273L |274| . The shell-model and nuclear structure calculations are crucial to compare theory 
to experi ments. The plethora of nuclear structure tech niques in volve, among others, the int eract ing shell model 
[275l l276| , the quasiparticle random-phase approximation |277L l278j| , the interacting boson model |279l | and the energy 
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density functional method (Mum. 


V. NUCLEAR ASTROPHYSICS 
1 . The Repulsive Coulomb 

Although stars are very hot, the kinetic energy of particles in their plasma are very small compared to the Coulomb 
repulsion energies, or the Coulomb barriers, in nuclei. For example, in the Sun’s core the temperature translates to 
about 10 keV of relative energy of the nuclei. It is very difficult to measure reactions at such energies because of the 
smallness of the cross sections due to the Coulomb repulsion. In general, i.e., either for the Sun or for other stellar 
scenarios, many of the reactions of interest are not known at the level of accuracy required by stellar modelers, or 
they are simply not measurable directly at all. Many reactions also involve short-lived nuclei which are difficult, if 
not impossible, to manipulate. 

One option for experimental nuclear astrophysicists is to use, e.g., stable or radioactive beams and study reactions 
with inverse kinematics techniques. Reactions with unstable nuclei use a stable target (e.g., a proton gas target) and 
by inference one can access information on cross sections of astrophysical interest. This information relies heavily on 
reaction theory |l35j |. We cite a few phenomena (only a few, indeed) of current interest using these techniques. 


A. Collective Resonances, Neutron Skins and Neutron Stars 


At low energies, of a few MeV/nuclei, heavy nuclei are not able to surpass the Coulomb barrier and nuclei basically 
interact only via the Coulomb force. At high energies, above the Coulomb barrier, experimental techniques have been 
developed to extract Coulomb excitation events from other sort of concurring reactions. A plethora of processes are 
studied using this technique which has the advantage of using a well-known interaction and therefore the reaction 
mechanism is very much under control. A revie w of the applications of Coulomb excitation to several reactions of 
astrophysical interest is found in Refs. |282l - l284 |. One of the processes of interest for astrophysics is the excitation 
of collective resonances above the particle emission threshold. The most notable of these are the so-called giant 
resonances. They exist in all sorts of modes. There are monopole, dipole, quadrupole and octupole vibrations, related 
to the character of nuclear shape vibration. They can also be of isoscalar or isovector character, due to their isospin 
vibration type (T = 0 ot T =1). There are also electric and magnetic collective excitations. 

Collective vibrations in nuclei are intimately related to the compressibility of nuclear matter, or equation of state 
of nuclear matter. This equation give the density dependence of the energy density in nuclear matter as 


E{p,a) = E{p,Q) + S{a)a^ + 0{a^), where a=— —(24) 

P 

with p = Pp + Pn equal to the nuclear density. The energy e{p) = E(p, 0) is the energy for symmetric nuclear matter. 
From it one can calculate quantities of interest for astrophysics, such as the pressure, P{p) = p^de{p)/dp, and the 
matter incompressibility, K — 9dP{p)/dp. These quantities, together with the entropy density, are key to determine 
the relation between mass and radii of neutrons stars. 

The EOS of nuclear matter can only be inferred by a intimate combination of theo ry and experiment because 
nothing on Earth can resemble the matter at extreme densities inside a neutron star [285j| . The main part of the EOS, 
i.e., the dependence of nuclear pressure on nuclear density can only be inferred in a small range of interest for neutron 
star modelers. These are often strongly guided by theoretical mean field calculations. Traditionally, theorists are able 
to calculate the nuclear matter energy incompressibility modulus, K, using energy density functional theories, which 
start from the nucleon-nucleon interaction in the nuclear medium and calculates E{p, a) using a microscopic theory 
such as the Hartree-Eock, Relativistic Mean Eield Theories, Random Phase Approximation, or variations of these 
theories [rMTll FT^. [28g. 

It is well-known that although some effective mean field interactions can reproduce many nuclear properties quite 
well, their pred ictio n strongly deviate from each other as one departs from the nuclear density saturation {p = po ~ 0.17 
nucleons/fm^ 1 [287| . As an example of such microscopic calculations, by comparing Hartree-Fock-Bogoliubov results 
for the excitation of giant monopole resonances in nuclei one can infer what sort of density functionals are most 
appropriate to reproduce the experimental data. Giant monopole resonances are the most direct compressible mode 
of a nucleus and and clearly related to the compressibility modulus. Inspecting t he m ost accepted forms ot the 
microscopic interactions has shown that a compressibility modulus, K ~ 200 MeV |l23j can be obtained, although 
values of AT ~ 230 MeV are also acceptable by studies of additional mean-field calculations for the excitation and 
decay of other giant resonance modes. 
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Relativistic heavy ion collisions also induce a matter compressibility during a short-time interval, with consequences 
such as side-flow of nuclear matter which can be measured and compared to micr oscopic c alculations usually done 
with transport theories such as the Boltzmann-Uehling-Uhlenbeck (BUU) equation [2421 l288l |. 


dl 

dt 


+ + VpC/^ • Vr/ - VrU • Vr/ = J Sp2 J a^N i^) |Vl - V 2 | 

X [1 - /i] [1 - /2] - / 1/2 [1 - f[] [1 - m ■ 


(25) 


Here, the number of particles in the volume element cPrcPp at time t, is given by / (r, p, t) (PrcPp, in terms of the 
distribution function / (r, p,t). The mean-field, U (r,p,t) accounts for the effect of each particle interacting with all 
others, with forces on each particle given by —V^U (r, p,t) and — Vpl7 (r, p,t) . The distribution function changes 
due to nucleons leaving (or entering) the volume d^rd^p, accounted for by the collision term on the right side, where 
aNN is the nucleon-nucleon scattering cross section, vi and V 2 are the velocities of two colliding nucleons, and their 
scattering angle. The factors (1 — /) account for Pauli blocking of final occupied states. 

In further constrain the EOS of nuclear matter, the function S in Eq. (1241) . called by the symmetry energy is needed. 
It is usually expanded around a: = 0 where x = (p — po)/3po, yielding 


1 ^ 
2 da^ 


a=0 


— J Lx ~h sym^ ' 


(26) 


where J is the bulk symmetry energy, L determines the slope, and K^ym the curvature of the symmetry energy 
at saturation, p = po (~ 0.17 nucleons/fm^). The use of heavy ion central collisions and the interpretation of 
experimental data using transport equations has been proven to be a good method to study the effects of symmetry 
energy in nuclei and in nuclear matter (see, e.g., Refs. 24,11 - 1245 11. 

Collective vibrations of nuclei are also a good tool to investigate the symmetry energy. Some features of neutron 
rich nuclei are also though to be very sensitive to the th e symmet ry energy, such as their neutron skin Ar = rp — r„. 


where rp(„) is the proton (neutron) matter distribution [2891 l290l| . C ollective vibrations in neutron-rich nuclei at low 
energies are usually so-called Pigmy Resonances or Ikeda resonances [IM[TM[29i[j^ . According to equation (14) 
of Ref. |294| . the energy of pigmy resonances, Epr, are elated to the neutron skin, Ar, and the neutron excess, Nexc, 
by means of 


EpH 


J. 




V 2Ari?TOjvA^e: 


1/2 


(27) 


where R is the nuclear radius and Sn the neutron separation energy. Inserting typical values for 5'„, R, Ng^c and Ar 
for medium mass neutron-rich nuclei, one obtains Epp = 5 — 8 MeV. 

Study of pi gmy reso nances helps constrain the symmetry energy dependence of the EOS of nu clear matter in 
neutron stars [2421129,5| . They are also related to the neutron skins in nuclei, as shown in Eq. dill). [28i,[29i,[293. 
Pigmy resonances are a possible energy sink during supernovae explosions. Their existence can lead to noticeable 
chan ges in the abundance of heavy elements by means of rapid neutron capture reactions, or r-processes for short 
[298l| . In the process a neutron capture occurs on a very short time-scale, ^ 0.01 — 10 s. This is often much faster than 
accompanying /3-decay processes occurring between the neutron captures. The r-process path therefore approaches 
the neutron drip-line as nuclei get more neutron-rich. A large number of isotopes in the range 70 < A < 209 are 
produced in this way. The stellar site of r-processes are not well known. It might occur during supernovae explosions, 
or during neutron star mergers processes. 

The nuclear dipole polarizability ap has recently been suggested [2991 - 1301 1 as an alternative observable constraining 
both neutron skin and symmetry energy. The polarizability is related to the photoabsorption cross section 


ap = 


he f aabs(uj) 


duj, 


(28) 


where uj is the photon energy. Because of the inverse energy weighting, the low energy response sensitive to the 
symmetry energy is probed. Recently, experiment s carried out in Japan and led by Peter von Neumann-Cosel allowed 
to infer the dipole polar izability of several nuclei [3021 - 1304 1. Early work by Takashi Nakamura and Thomas Aumann 
[l89l I 191 LI 192 L l29ll . |292|| have also unravelled the electromagnetic response of pigmy resonances with gain of insight 
into the symmetry energy. 
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B. Nuclear Reactions in Stars 


1. Reaction Rates and Radiative Capture 


In stars most nuclear reactions occur in the form of binary collisions. One of the exceptions is the triple-alpha 
reaction mentioned in the introduction. For nuclei j and k in an astrophysical plasma, the reaction rate in binary 
collisions is given by = (av) rijUk, where the reaction rate {av) is the average of the product of the cross section 
and the relative velocities, av, over the temperature distribution. More specifically, 


/ \ VljTlk 

r,- 1 . = (av) ■, — 

\ h,k I g 


'jk 


(29) 


where 




TTljj^lT 


1/2 


Ea{E) exp dE. 


(30) 


Here rrijk denotes the reduced mass of the target-projectile system. The factor 1 -|- Sjk in Eq. (1^ accounts for the 
case of identical particles, i = k. The rate should be proportional to the number of pairs of interacting particles in 
the volume. If the particles are distinct, that is just oc njUk- But if the particles are identical, the sum over distinct 
pairs is oc rijUk/^. 

In Eq. (1501) the velocity, or kinetic energy E, distribution of the nuclei in the plasma at temperature T is assumed to 
be Maxwellian. It has been shown in many instances that even a small depart from the Maxwell-Boltzmann velocity 
distribution would modify dramatically the outcome of the stellar evolution. For example, recently it has been shown 
that alternative statistics would ultimately have to be extremely close to Boltzmann-G ibbs statistics to reproduce the 
observed relative abundance of elements evolved during the big bang nucleosynthesis |305| . 

The direct measurement of the reaction cross sections a for nuclear astrophysics is difficult because it either involves 
neutron-induced reactions such as in the r-process, or charged particle reactions which are strongly inhibited by the 
Coulomb repulsion, or Coulomb barrier. Many of such reactions also involve nuclei far from the stability or neutron- 
rich nuclei. These nuclei are not usable as targets in laboratory experiments. Thus, most experiments with nuclei far 
from the stability requires the use of radioactive beam facilities. They are also done at much higher energies than the 
typical nuclear relative motion energy in the stars. Therefore, many indirect methods have been devised to extract 
indirectly the cross section for nuclear reactions in stars. 

One of such indirect methods is known as the Coulomb dissociation method. This method is useful when a projectile 
is loosely-bound and can be dissociated (a —)> 6 + c) in the Coulo mb field o f a target with a large charge Z. The 
Coulomb breakup cross section for a -I- H & -|- c -|- H is given by [loillQl, 


dag 

dEdO. 


^ dEdn ^ b+ciE), 

ILL 


(31) 


where H = ^(electric) or M(magnetic), T = 1, 2, 3, • • • is the multipolarity, and ^ b+ciE) is the the photonuclear 
cross section for photon energy E. dN^^/dEdfl are the equivalent photon numbers, depending on E, and the s cattering 
angle = {9, cj)) for a -I- H —5- b + c + A. The equivalent photon numbers are obtained theoretically [306j | for each 
multipolarity HL. Time reversal allows one to relate the radiative capture c ross section 6 -|- c —> a + 7 to the 
photo-breakup cross section ^ f,+c(w). This method was proposed in Ref. [307j and has been applied to several 
reactions of interest for astrophysics. For example, the reacti on ^Be(p, 7 )®B relevant for the production of high-energy 
solar neutrinos has been studied with this method [3n3 - l3Tnl| . 

For many systems studied via the Coulomb dissociation method, the contribution of the nucle ar b reakup is relevant 
and cannot be ignored. Such contributions has been examined by several authors (see, e.g. [31l| l. In t he case of 
loosely-bound systems such as ®B it has been shown that multiple-step, or higher-order effects, are important |312l - l314 
due to continuum-continuum transitions. The fragments can make excursions in the continuum and iterate with the 
Coulomb field of the target many times before they become asymptotically free. Multistep processes in the continuum 
led to the development of new techniques, such as th e continuum-discretized coupled-channels (CDCC) method, e.g., 
as applied to a loosely- bound nu cleus in Ref. |312| |. The effects of high multipolarities (such as E2 contributions) 
have also been studied [3I5l - l317 | for t he reaction ^Be(p, 7 )®B). It has also been shown that the influence of giant 
resonance states can be neglected [3I8j| . The challenge is to develop reliable nuclear models to obtain the scattering 
states including resonances. 















17 


2. Single-Particle, Ab-Initio Methods 


Taking radiative capture reactions as an example, the cross sections for direct capture are given by (313 

^Elj,^\{lcJc\\OnL\\hjb)\\ (32) 

where a proportionality factor involving phase-space quantities is omitted, Oul is the electromagnetic operator, and 
(lcjc\\OuL\\hjb) is a multipole matrix element involving bound (6) and continuum (c) wavefunctons. For electric 
multipole transitions and in the long-wavelength approximation, OnL = 

nOO 

{Icjc \\Oel\\ hjb) OC / dr Ub{r)uc(r) , (33) 

Jo 

where Ui are radial single-particle wavefunctions. If the wavefunctions are obtained in a many-body model, including 
anti-symmetrization and correlations, the equations can get more involved. If a relation with the single-particle model 
can be made, then the total direct capture cross section is obtained by adding all multipolarities and final spins of 
the bound state {E = Enx), 


L,Jb 


(34) 


where (S'F)j,, are spectroscopic factors. The spectroscopic factors are adjusted to a shell model calculation for the 
single-particle occupancy amplitudes. This kind of calculations have been rou tinely used in the literature for the 
^Be(p, 7 )®B) and other radiative capture reactions (see, e.g., Refs [3171 l3I9l - l324 |l. 

Radiative c apture re actions using EFT techniques, such as those mentioned in an earlier section have been reported 
in e.g. Refs. |325l - l327l | . The advantage of this approach is the predictive power of EFT which allows for increasing 
the precision of calculations as more terms in the effective theory expansion is included (the dots). One of the main 
problems of this method is the inclusion of many-body intrinsic features such as anti-symmetrization and medium 
modification of the interactions. 

The Resonating Group Method (ROM) or the Generator Goordinate Method (GCM) are able to include anti- 
symmetrization of bou nd a nd continuum states on equal foot. These are based on a set of coupled integro-differential 
equations of the form [328l| 


J d?r' [H^l^,{v,v') - EN^^,{v,v')\ ^^/(r') = 0, 


(35) 


where E[{N){y,y') = ('I'^(a, r)|iJ(I)|'I'B(a', r')). Here, E[ is the Hamiltonian for the system of two nuclei (A and B) 
with the energy E, a,b is the wavefunction of nucleus A (and B), and gaiO is & function to be found by numerical 
solution of Eq. (1351) . This function describes the relative motion of A and B in channel a with full antisymmetrization 
between nucleons of A and B. Such calculations have been carried out successfully for many reactions of astrophysical 
interest (see, e.g.. Ref. [13). 

To the present date, the RGM is the most practical tool to describe nuclear scatt ering and fusion of light nuclei 
at low energies with the inclusion of anti-symmetrization. It was introduced in Ref. [330| and by solving it one can 
calculate cross sections for reactions involving cluster-like nuclei such as ^^G(a,7)^®0 (33ll |. This reaction is relevant 
for cosmology. It is thought that if the cross section for this reaction would be twice the adopted value, a 25 solar 
masses star will not produce ^°Ne since carbon burning would cease. An oxygen rich star is more li kely to collapse 
into a black hole while carbon rich progenitor stars is more likely to leave behind a neutron star |332|. Can you 
imagine how important this nuclear reaction is for Cosmology? This reaction determines the density of black holes 
and neutron stars scattered through the Universe. 

Other techniques, beyond the single-particle and RGM models have been pursued. The ultimate goal is to calculate 
the cross sections from an ab-initio-type of calculations, i.e. a calculation that starts from the bare nucleon-nucleon 
interaction an d evo l ves t o solve the many-body nuclear problem, including the continuum. Eirst attempts have been 
done in Refs. [33311334| with ab-initio bound state wavefunctions calculated with the Quantum Monte-Carlo Shell 
Model or the No-Core Shell Model. In both cases the continuum wavefunction was derived by assuming a single¬ 
particle cluster mo del for th e scattering states. Later, the RGM has been combined with NSCM to calculate radiative 
capture reactions |335l . l33(| . Another example is a calculation of the tr ansfe r reaction ^H(d,n)‘*He cross section to 
explain the experimental data obtained at the National Ignition Facility |337j . 
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3. Transfer, Trojan Horses, ANCs, and Surrogates 

Although some direct measurements are possible, calculating transfer reactions such as ®Li(d,a)^He is relatively 
difficult. In general, the measurements also suffer from the same kind of difficulties we mentioned earlier. Therefore, 
some indirect methods have been developed. Among these we cite the Trojan Horse Method (THM). Because of 
off-shell effects a nucleus x can be carried inside a trojan nucleus and brought to react at low energ ies with a heavy 
target. The heavier and faster trojan nucleus allows x to overcome the Coulomb barrier [338l l339l | . The basic idea 
is to extract the cross section in the low-energy region of a two-body reaction with significant astrophysical impact, 
a + X —> c + C, from a suitable three-body reaction, a + b — > s -|- c -I- C, at high energies. An example is the study 
of the S-factor for the ^H(d,p)^H reaction by using ^He breakup in the ^H(^He,pt)H reaction. The analysis of THM 
data is done by applying the well-known theoretical formalism of the quasi-free (QF) process. One of the clusters in 
the incoming projectile is the participant, x (the deuteron in the given example), while the spectator s, (^He or p) 
are only weakly involved in the process. These assumptions allows one to relate the cross section to the momentum 
distribution of the inter-cluster (x — s) motion inside b 

Trojan Horse (TH) reactions allow for the use of higher bombarding energies, Ea, high enough to overcome the 
Coulomb barrier of a -I- a: in the entrance channel of the reaction. Then the effect of the Coulomb barrier and electron 
screening effects become small. The triple differential cross section in the center of mass of the TH reaction can be 
written as 


d^c 


dEcdTlcd^c 


= Kf |$(p.x)| 


h 


daj, 


dflr 


HOES 


(36) 


where k is the orbital angular mome ntum of particles s and x and is a function of relative momentum and 
kinetic energy in the entrance channel |34Cll| . In this equation, {da /is the half-off-energy-shell (HOES) 
differential cross section for the two-body reaction at the center-of-mass energy given by Ecm = Ec-c — Q 2 b, where 
Q 2 b is the two-body Q-value of the binary process and Ec-c is the relative energy between the outgoing particles, Kp 
is a kinematical factor, and ^{psx) is the Fourier transform of the radial wave function y {r) for th e x — s intercluster 
motion. The success of the THM relies on the QF kinematics (equivalent to d)(ps 2 ;) ^ 0 |34ll - l345j | . 

One h as sh own that the THM obtains is the same cross section energy dependence as those obtained with direct 
methods |346|| . The THM data are renormalized to those obtained with direct methods at high energies and the data 
at low energies are obtained as a bonus. The THM has been applied to a large number of reactions of inter est for 
the production of light elements in numerous stellar environments to reactions of interest for the BBN |347j| . One 
of the many advantages of usin g transfer reaction techniques over direct measurements is to avoid the treatment of 
the electron screening problem 348l - l35(3j| . This is just because the reaction is peripheral, probing the tails of the 
bound-state wave functions. In fact, most astrophysical nuclear fusion processes reed only this information for an 
accurate theoretical calculation. 

The Asymptotic N ormalization Coefficients (ANC) method, introduced by Akram Mukhamedzhanov and collabo¬ 
rators |349l . l35lL l352j| , assumes that the amplitude for the radiative capture cross section b + x —> a -I- 7 is calculated 
according to 


M = (ibxi'^b^) O(rbx) V'l’^^(rbx) 


(37) 


where 


^bx — \^a{£,bj fxj Tbx) 4>x{f,x)4’b{£,b) 


(38) 


is the overlap function, i.e., an integration over the internal coordinates and fx, of b and x. At low energies, is 

dominated by contributions of large rbx- Thus, the matrix element M is dominated by the asymptotic value of 


^bx ^bx ^-r)a,l/2(2Kba:?'ha;)) 
'^bx 


(39) 


where Wa,n is the Whittaker function and is the ANC. The ANC can be related to single-particle properties by 
writing it as a product of the single particle (s.p.) spectroscopic factor and a normalization constant which depends 
on the details of the s.p. wave function in the interior part of the nucleus. If the reaction occurs at very low energies, 
then the radiative capture cross section only probes the tails of the bound-state wave functions. Therefore, is the 
only unknown part of the bound state wavefunction needed to calculate the direct capture cross section. The ANCs 
can be obtained by experimental analysis of elastic scattering between nuclei by extrapolation of the scattering phase 
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shifts data to the bound state pole in the energy plane. They can also be accessed in peripheral transfer reactions 
whose amplitudes contain the same overlap function as the amplitude of the corresponding astrophysical radiative 
capture cross section. 

As an example, consider the proton transfer reaction A{d, a)B, where d = a + p, B = A + p. Using the asymptotic 
form of the overlap integral, the cross section is given by 


da 


E 

JBjd 




1 — 

(N 

s 

1_ 

L \ 


P'^ 


(40) 


where ct is a reduced cross section. Pap {Pap) is the asymptotic normalization of the shell-model bound-state proton 
wave functions in nucleus d{B), which are related to the corresponding ANC’s of the overlap function by means 
of [CapY = SapPapi where S^p is the spectroscopic factor. If the reaction A(d,a)B is peripheral, the bound-state 
wavefunctions used to calculate d are also approximated by their asymptotic form, leading to ct oc P\pPap- Hence 

%=Y^{CApf{CtpfRBd, where Rsd = 

““ PApPap 


is independent of P\p and Pap- Thus for surface-dominated reactions, the cross section is obtained in terms of the 
product of the square of the ANC’s of the initial and the final nuclei (C'^p)^(C'fp)^ rather than in terms of spectroscopic 
factors. 

Most of the (n, 7 ) reactions of astrophysical interest will never be measured directly. Transfer reaction methods as 
the ones described a bove are form a kind of surrogate reactions in which the neutron is carried along by a nucleus to 
react with a target |353l l354l |. An example of surrogate reactions are (d,p) reactions, which exploits the simplicity 
of the deuteron structure. Evidently, the neutron in a nuclear environment carries angular momentum and other 
quantum numbers which are different than those of free neutrons. Unless the angular momentum plays a minimal role 
in the reaction, this brings difficulties in extracting the desired neutron-induced reaction from the surrogate equivalent. 
Theoretically, this is the same a s the implication that the Hauser-Feshbach formalism for ( 0 , 7 ) reactions agrees with 
the Ewing-Weisskopf formalism |355|| . To elaborate on top of our discussion of the Hauser-Feshbach formalism leading 
to Eq. (11311 . assume that (TAr(c) is the cross section for the formation of a compound nucleus in the entrance channel 
c. Using the time-reversal theorem, the cross-section dec' can be related to the cross-section for the time-reversed 
process c' c, and one gets 


Ccc' {Ec')dEc' = acN (c) 


(2/c' -I- 1) pc'Ec' acN {c')u}{Uc')dUc' 

_ „ rmax 5 

fo ° (2-fc + 1) dcEc acN (c) oj{Uc)dUc 


(42) 


where R is the angular momentum and pc is the reduced mass in channel c. Fragments emitted with energy in the 
range Ec' to Ec' +dEc' leave the residual nucleus with energy in the range Uc' to Uc'+dUc' where Uc' = Ecn — Bc'—Ec' 
and Ecn and B^' are the compound nucleu s ene rgy and the binding energy of the fragment in the compound nucleus. 
Eq. (14^ is the Weisskopf-Ewing formula |144| . The level density is approximately given by io{U) oc exp{U/T), so 
that the emitted fragments in the Weisskopf-Ewing theory follow a Maxwell evaporation spectrum. 

The Weisskopf-Ewing theory does not explicitly consider the cons ervation of angular momentum and parity J and 
parity tt. This is described by the Hauser-Feshbach expression |l43l| 


TT (2J-fl) Es.j^*(c)Es',i''H/'(c') 

P 4" (2*c + 1) (2/c + 1) EeEs,lTl{c) 


(43) 


where T/ are barrier transmission probabilities that depend on spin and parities. If the (d,p) reaction is well reproduced 
by the E wing -Weisskopf theory, this is a hint that the reaction is a good surrogate to obtain (n, 7 ) and (n,fission) cross 
sections |356l |. This means that the surrogate reaction does not depend on spin and parities and that the (d,p) and 
(n, 7 ) and (n,fission) reactions populate the same states. 

Surrogate reactions involving charged particles, such as (d,p) reactions are complicated because they require a 
correct treatment of the interact ion b etween its constituents to all orders. A popular tool to treat these systems is the 
Alt-Grassberger-Sandhas (AGS) [357j| formalism. They arise from the Faddeev formalism treating the system {a-\-b)-\-A, 
and its rearrangements a-\-{h-\- A) and h-\- {a-\- A). In each of these channels the three-particle scattering is described 
in terms of the transition ope rators Tpa, where ol{P) corresponds to a set of channel permutation combination. They 
obey the AGS equations |357l | that are a system of coupled integral equations given by 

3 

7>a(z) = -f ^ 5l}vtu{z)Go{z)T^a{z) , 


(44) 
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where = 1 — Sp^, is the anti-delta-Kronecker symbol, Gq = {E + i^+ — Ho)~^ is the resolvent of the three-free 
particle c.m. energy for the Hamiltonian Hq. The two-body transition operator for each interacting pair with 
inter-potential v^, is obtained from the Livymann-Schw inqe r equation + v^Got^. The solution of these set of 

equations require some approximations (see, e.g.. Ref. 358l |') and one of the main problem is the accurate treatment 
of the Coulomb interaction. Some methods have been developed to handle the Coulomb interaction by means of 
screened potential, but i t onl y works for relatively light nuclei because the solutions oscillate wildly at the asymptotic 
region for large charges [35^. 

During supernovae core collapse, temperatures and densities are so high that electrons are captured in nuclei 
reducing the electron fraction in the medium, and driving the nuclear composition to more neutron rich and heavier 
nuclei. As a consequence, nuclei with TV > 40 dominate matter composition for densities larger than a few 10^° g 
cm“^. Simulations indicate that electron capture stops at such densities and the capture becomes entirely due to free 
protons. In order to ensure that such simulations are correct, one needs to study electron capture in a large number 
of nuclei. These studies cannot be carried out in a direct manner in the laboratory. The main object of interest 
are Gamow-Teller matrix elements, B{GT), which cannot be extracted from /3-decay experiments The method of 
charge exchange reactions was developed to access the B(GT) matrix elements using (p,n), (^He,t) and other charge 
exchange reactions at energies around 100 MeV/nucleon [360j. In such reactions one obtains the nuclear response 
by spin and isospin operators also involved in neutrino-s catte ring reactions occurring in stellar environments, useful 
to understand the neutrino driven explosion mechanism [361| . This approach relies on the similarity in spin-isospin 
space of charge-exchange reactions, electron capture, neutrino s catte ring, and (3 decay mechanisms. For example, the 
cross section <j{p,n) at small momentum transfer q is given by [362| |. 


da 

dq 


{q = 0) = KND\J.r\‘^Bia), 


(45) 


where AT is a kinematical factor. No is a distortion factor (accounting for initial and final state interactions), J^t is the 
Fourier transform of the effective nucleon-nucleon interaction, and B{a = F, GT) is the reduced transition probability 




.(±)i 


for non-spin-flip, B{F) = {2Ji + 1) 
transitions. 

The charge-exchange matrix element is given by [363 


f)P, and spin-flip, B{GT) = {2Ji + 1) ^|(/|| ||f)| 


Mexchiq) = 








(46) 


(if) 

where q is the momentum transfer, '1')^ ^ are the intrinsic wavefunctions of nuclei a and b for the initial and final 
states, Va^b are the nucleon coordinates within a and b, and Vexch is the part of the nucleon-nucleon interaction 
containing spin and isospin operators. For forward scattering and low-momentum transfers, q ~ 0, and the matrix 
element (l46l) becomes 


Mexchiq - 0) - u® Jq - 0) Ma(F, GT) Mb{F, GT) , (47) 

where idoes not depend on spin-isospin, and Mexch{F,GT) = or cr)T|are Fermi or Gamow- 

Teller (GT) matrix elements for the nuclear transition. A similar result can be obtained using eikonal scattering 
waves instead of plane waves. This justifies the use of Taddeucci’s formula, as long as one remains withi n th e 
validity of the low-momentum transfer assumption in high-energy collisions, as discussed thoroughly in Ref. [363|| . 
The prop ortio nality between cross section and Gamow-Teller strength was shown to apply to a wide range of mass 
numbers |364j| . 


4- Hadrons, Quark Soups, and Strong E/M Fields 

It is predicted that during the Big Bang, a phase transition occurred from a soup of quarks and gluons to a hadronic 
phase. This soup was present in an environment of high temperature and densities. A phase transition from hadrons 
to quarks and gluons, although in a cold environment, might also occur in the core of neutron stars. The only possible 
way to study phase transitions to sub-nucleon degrees of freedom are by means of ultra-relativistic heavy ion collisions 
(URHIG), routinely carried out at the RHIC accelerator in Brookhaven and at the LHC at CERN to induce very high 
temperatures, densities, and pressures in central collisions. Under such conditions, it is expected that the nucleons 
are excited into baryonic resonance states, along with particle production crea ting a hadronic resonance matter and 
a possible transition to a free soup of quarks and gluons during a short time [sj, l365[ | . QCD induced phenomena, such 
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as quark confinement and spontaneous chiral symmetry breaking are caused by infrared physics based on asymptotic 
freedom, but it is not fully understood how they are related with each other. LQCD simulations predict that these 
two phen omena are apparently turned off at about the same time at a given temperature at the onset of the QGP 


By studying URHICs one extracts of the energy density e, pressure P , and entropy density s of the soup as 
a function of the temperature T and the baryochemical potential ps- During a phase transition to a quark-gluon 
plasma (QGP) a rapid rise in the effective number of degrees of freedom, expressed by e/T^, or s/T^, is observed. 
The variables T, s, and e, are correlated with the average transverse momentum {pt), the hadron rapidity density 
dN/dy, and the transverse energy density dEx/dy. The transverse energy produced in the interaction is given by 
Et = j where Ei and 9i are the kinetic energies of the fragments and their emission angles. Photons 

and leptons produced in URHIG provide information on the various stages of the interaction without modification by 
final state interactions. The widths and masses of the mesons p, uj, and (j) appearing in the lepton pair invariant mass 
spectrum are expected to be sensitive to medium effects. Jparticles are abundantly produced in a quark-gluon 
plasma and predicted to be suppressed due to the Debye screening of the cc pairs. Loosely-bound states of the cc 
system, such as d/' and Xc , are more easily dissociated and are even more suppressed than the J/d>. In a QGP the 
strange quark content is rapidly saturated by ss production in gluon-gluon interactions. This leads to an enhancement 
in the production multi-strange baryons and strange antibaryons when a QGP is formed. 

The energy loss of a quark is related to the color-dielectric polarizability of the medium, in analogy with the 
energy loss by electromagnetic probes. Due to Bremsstrahlung, radiation is a very efficient energy loss mechanism 
for relativistic particles. But it is strongly suppressed in a dense medium by the Landau-Pomeranchuk effect Q. By 
analogy, the stopping power of particles within a QGP is predicted to be higher than that in the hadronic matter, 
leading to an attenuation the emission of gluon jet pairs from parton-parton collisions in the direction opposite to the 
trigger jet, known as jet-quenching and visible in measurements of high-pr spectra. A quark or gluon jet propagating 
through a dense m edium will not only loose energy but will also be deflected, leading to an azimuthal asymmetry of 

the jets [inl-izi- 


In order to understand the QGP phase transitions, the boundary conditions in the collisions have to be understood 
properly. For example, the parton distribution functions (PDF) which are basically medium modified gluon distribu¬ 
tions in nuclei have to be stu died carefully. The strong electromagnetic field generated by one of the colliding ions has 
been proposed in Ref. [376j| t o be used as a tool to study PDFs at the LHG/GERN. The physics of ultra-peripheral 
collisions is reviewed in Refs. The sensitivity of the different modern PDFs to t he mediu m corrections 

is clearly visible in the rapidity spectrum of J/d^ a nd T particles produced via this process [Izi-lilll. The power 
of experiments to discern among the several PDFs l379l| th ose that best reproduce the experimental data has been 
evidenced in experiments at the LHC/CERN |l23L l38 J 383 1. This allows for a better discrimination of the pre-initial 

PDE conditions for interest in the formation of the QGP in URHICs. _ 

The experimental discovery of a QGP was announced in different occasions in the last few decades [3841 - 1388 1. 
The numerous press releases with announcements of the QGP discovery reflects the extreme difficulty to produce 
and confirm the existence of such an ephemerally state of nuclear matter, lasting an extremely short time with 
temperatures as high as 10^^ K. Of present interest is also the glasma state of matter occurring when each of the two 
nuclei in URHIG form a color glass condensate, related to the PDEs, preceding the QGP. In such a state the dressed 
parto ns are con densed into a glassy state, prior to the transition between the confined state within hadrons and QGP 

UliMIl- 

After the experimental discovery of the Higgs boson |392l . l393| , the CERN experimentalists are searching for new 


state 


phenomena beyond the stand ard model of particle physics, among others there are dedicated searches for super- 
symmetric particles [3941 - 1398 1 and extra-dimensions [399l - l403| . But there is also a very strong community at CERN 
studying nuclear physics, not only with regard to the formation and decay of the QGP, but also for other purposes. At 
present, exciting experiments at m uch lower energies are being carried out in nuclear physics at the CERN/ISOLDE 
radioactive ion beam facility [404| and many other sort of breakthrough experiments in nuclear physics make the 
headlines now and then . A few exam ples are th e early production of anti-hydrogen, the first time and anti-atom was 
produced and detected Hoi (see also [iofil - lin^ . or the recent precision measurements of the mass difference between 
light nuclei and anti-nuclei [ 403 . 


VI. DREAMS COME TRUE 

One of the dreams of human technology is to reproduce on Earth the same fusion processes occurring in the Sun, 
namely, the immense energy release through fusion of light nuclei. We have developed Tokamaks, marvelous machines 
that conhne the hot plasma necessary to trigger nuclear fusion. However, we have never been able to develop a machine 
that produces more fusion energy than it is employed to run it. It is a pity, as fusion energy is free from many of 
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the problems related to nuclear reactors, most notably the immense nuclear radiation problem. Large multinational 
projects seem to be the only way to couple with such incredibly expensive tasks. We are right now building a fantastic 
Tokamak machine in France, with funds from several countries flowing in. The project is very costly, possibly reaching 
20 billio n Eu ros or more by the end of its completion. As of today, the machine is predicted to obtain the first plasma 
by 2030 |410l| . Some of us will not be here to witness it. 

Unfortunately, ITER (pronounced “eater”), an acronym for International Thermonuclear Experimental Reactor, 
is running into trouble, as countries funding the project do not deliver their sub-projects on time, delaying the 

a rammed deadlines. A similar case, at much smaller funding scale, is going on with the FAIR facility in Germany 
. It seems that the best approach for expensive multinational projects is to allow a single country to run most 
(if not all) of the project steps, much like the CERN model has proven to function with success. From the science 
point of view, ITER does not plan much of nuclear science research. Basically the project is related to material 
science and the studies of what radiation can do to the material confining the plasma. A new ITER possibly built 
many decades after (we will all be dead), might then reach the so much expected fusion energy gain. This is sort of 
disappointing. But it might also be a fallacy of the way we think about such projects. Who would predict 3 decades 
ago that studies of anti-hydrogen might be of interest for the CERN laboratory?, or that tests of ab-initio nuclear 
structure calculations might be of interest for the National Ignition Facility (NIF)? In the end, nuclear physicists will 
certainly find their niche in the ITER project, just because nuclear physics research is resilient and necessary. 

Many other dreams are out there, some of them perfectly doable within present fundi ng c apabilities, or already 
being accomplished. Among a few, we cite heroic efforts by Toshimi Suda leading SCRIT (41 ij . a novel experimental 
technique to trap short-lived nuclei almost at rest, a completed electron scattering facility inside RIKEN/RIBF 
employing a recycled electron storage ring. First experiments are underway. Also worth mentioning are the efforts 
of Raju Venugopalan, Abhay Deshpande, and collaborators to promote the idea of an Electron-Ion Collider (EIC) in 
the U.S. using extensions of present facilities. Such a facility would probe the structure of nuclei down to distance 
scales as small as 10“^ fm allowi ng th e study o f mo mentum, spatial, spin and orbital distributions of gluons and sea 
quarks in light and heavy nuclei |412| (see also [4I3|| b 

Many similar ideas are flourishing in nuclear physics and many of them are found to be extremely useful for the 
progress of science. 


VII. CONCLUSIONS, EXCUSES AND COMPLAINTS 


In this review we have attempted to convey a realistic picture of the current status of nuclear physics research. 
The choice of the topics covered in the review reflects our own interest and contributions to the field. Clearly several 
topics were left out. The message which we hope to transmit is that the field is quite alive with future planning 
and developments are underway and will guarantee a thriving field of physics research in the next decades. Nuclear 
physics is a synergetic field in the sense that ideas, theories, and concepts find a natural application in other fields, 
such as cold atomic gases (Feshbach Resonances), atomic and molecular collisions (CDCC), and universal conductance 
fluctuations in mesoscopic systems such as open quantum dots and grapheme flakes (Quantum Chaotic Scattering 
and RMT). 

In the review we have also discussed the recent development in the reaction theory used to treat the collisions 
of weakly bound and halo nuclei. The fusion and incomplete fusion of these systems are discussed and the theory 
employed robustly analyzed. We have also discussed the relevance of nuclear physics to astrophysics and presented 
a rather detailed account of the type of reactions of importance to nucleosynthesis. The underlying QCD and its 
constraints on nuclear phenomena were also discussed. The quark-gluon plasma predicted to exist, albeit for a very 
short time, in the high LHC energy collisions of heavy ions is also briefly discussed. 

In summary, in this review we have tried to give a flavor of what Nuclear Physics means for science, and what are 
its most actual academic issues. There are more facets to Nuclear Physics than we can cover in such a short review 
and we had to leave out the discussion of a large number of other interesting research topics. 
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